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AN  EXTENSION  OF  A  KINETIC  THEORY  OF  POLYMER  CRYSTALLIZATION  THROUGH  THE 

EXCLUSION  OF  NEGATIVE  BARRIERS 


Jerry  I.  Scheinbeim,  Louis  Petrone,  Brian  A.  Newman 
Department  of  Mechanics  and  Materials  Science 
Rutgers  University 
Piscataway,  NJ  08855-0909 

ABSTRACT 

The  simplest  version  of  the  Lauritzen-Hoffman  (LH)  model  of  polymer 
crystallization  which  applies  to  infinitely  long  model  polymer  molecules 
crystallizing  on  an  existing  substrate  of  infinite  width,  is  re-examined.  The 
mathematical  expressions  for  the  model  free  energy  barriers  are  observed  to 
take  on  negative  values  at  high  supercooling.  Since  such  negative  barriers 
appear  to  be  physically  unrealizable  for  the  crystallization  process,  the  LH 
model  is  extended  by  imposing  a  mathematical  constraint  on  the  expressions  for 
the  barriers,  to  forbid  them  from  ever  being  negative.  The  extended  model 
contains  one  parameter  7  which  varies  from  zero  to  one  and  is  analogous  to  the 
parameter  ij>  of  the  LH  model.  For  all  values  of  7  less  than  one,  the  extended 
model  predicts  a  finite  lamellar  thickness  at  every  supercooling;  moreover, 
this  thickness,  at  large  undercooling,  decreases  monotonically  with  increasing 
undercooling  in  agreement  with  experiment,  but  in  marked  contrast  to  the  LH 
model  which  exhibits  the  well-known  Si  catastrophe.  The  relative 
insensitivity  of  the  calculated  lamellar  thicknesses  to  the  parameter  7 
supports  the  use  of  7  -  0  as  a  first  approximation  for  mathematical 
convenience  in  practice. 


1. 


I .  INTRODUCTION 

Recently,  the  crystallization  of  poly (vinyl idene)  fluoride  (PVF2)  in  the 
presence  of  high  electric  fields  has  been  studied  both  experimentally  and 
theoretically.  Of  the  four  well-known  crystalline  forms  a,  /?,  7,  and  S  (or 

II,  I,  III,  and  IV)  of  PVF2,  the  phase  with  the  largest  spontaneous 
polarization  and  potential  for  applications  is  the  ^-phase. 

Crystallization  of  PVF2  from  a  concentrated  solution  of  tricresyl 
phospate  in  the  presence  of  a  high  electric  field  was  observed^-  to  produce  /?- 
phase  crystals,  with  dipoles  oriented  in  the  field  direction,  during 

the  initial  stages  of  crystallization  followed  by  the  growth  of  unoriented  a- 
crystals  (non-polar)  as  crystallinity  increased  and  the  tricresyl  phosphate 
content  decreased  by  evaporation.  The  decrease  in  tricresyl  phosphate  content 
and  subsequent  crystal  growth  behavior  suggests  that  the  local  electric  field 
in  the  solution  region  changes.  Other  experimental  and  theoretical  ’  studies 
of  crystallization  of  PVF2  from  the  melt  in  the  presence  of  a  high  static 
electric  field  have  been  made,  and  were  found  to  give  7-phase  crystals  which 
however  did  not  show  crystal  orientation.  As  part  of  the  continuing  effort  to 
understand  the  structure-property  relationships  of  FVF2  and  because  of  its 
practical  importance,  our  ultimate  goal- -despite  the  complexity  of  the  system 
described- -  is  to  develop  a  theory  or  model  which  can  account  for  its 
crystallization  behavior  from  concentrated  solutions  in  the  presence  of  an 
electric  field. 

As  in  the  case  of  isothermal  crystallization  of  a  and  7  phase  from  the 

3 

melt  in  an  electric  field,  a  theory  of  isothermal  crystallization  of  a,  p, 
and  6  phase  from  concentrated  solution  in  an  electric  field,  would  be  based  on 
"classical"  and  "polymer"  theories  of  nucleation  and  growth  in  the  absence  of 
an  applied  field.  Most  importantly,  the  nucleation  barrier  or  activation  free 
energy  barrier  for  nucleation  would  certainly  be  different  in  the  presence  of 
the  field  than  in  its  absence;  and  this  barrier  has  been  seen  to  be  of 
fundamental  importance  in  the  theories  of  polymer  crystallization,  the 
simplest  of  which  is  the  LH  or  Lauritzen-Hoffman  theory.^’ ^  One  possibly 
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unrealistic  feature  which  seems  to  have  been  incorporated  into  the  LH  theory 
in  order  to  simplify  it,  is  that  the  nucleation  barrier  is  not  constrained  in 
the  theory  to  take  on  only  nonnegative  values.^  The  word  "barrier"  connotes  a 
positive  quantity,  and  furthermore,  the  LH  theory  is  based  on  transition  state 
theory  *n  which  the  barrier  corresponds  to  an  intermediate  configuration  or 

transition  state  of  the  system  which  is  at  a  free  energy  maximum  relative  to 

9 

some  initial  and  final  state  of  the  system.  Moreover,  the  LH  theory 

exhibits,  in  contrast  with  experiment,  the  Si  catastrophe  wherein  the 

calculated  average  lamellar  thickness  2  suddenly  passes  through  a  minimum  and 

becomes  infinite  at  a  temperature,  T  ,  corresponding  to  a  moderately  large 

undercooling ;  and,  in  fact,  the  nucleation  barrier  in  this  theory  is  positive 

for  all  T  >  Tc,  is  zero  at  T  -  T  ,  and  is  negative  for  all  T  <  Tc  for  the 

special  case  which  Lauritzen  and  Hoffman  have  recently  considered. 

Therefore,  prior  to  developing  an  extension  of  the  LH  theory  which  would 

involve  ascertaining  the  effect  of  an  electric  field  on  the  nucleation 

barrier,  we  try  to  extend  the  LH  theory  to  larger  undercooling  by 

incorporating  into  it  the  assumption  that  free  energy  barriers  cannot  be 

negative.  Note  that,  unlike  in  the  LH  theory  of  polymer  crystallization, 

barriers  in  classical  nucleation  theory  are  never  negative;  however,  the 

classical  theory  does  not  explicitly  take  into  account  polymer  chain  folding, 

and  for  that  reason,  we  have  not  yet  considered  modifying  the  Marand  and  Stein 
o 

theory  of  crystallization  from  the  melt  to  treat  the  P^^/tricresyl  phosphate 
crystallizing  solution. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  II,  the 
LH  model  is  described.  The  kinetic  treatment  of  the  LH  model  is  given  in 
Section  III.  The  rate  constants  needed  for  this  treatment  are  determined  in 
Section  IV.  Next,  our  extension  of  the  LH  model  is  described  in  Section  V; 
the  conditions  which  determine  the  sign  of  A^,  the  free  energy  of  formation 
of  that  portion  of  a  model  polymer  molecule  which  crystallizes  first  on  an 
existing  crystal,  are  found  in  Section  VI.  A  summary  of  the  expressions  for 
the  barriers  in  cur  model  is  given  in  Section  VII  along  with  the  expressions 
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for  the  average  lamellar  thickness.  In  Section  VIII,  the  variable 
transformations  required  as  a  preliminary  to  numerical  integration  are 
introduced.  Results  and  discussion  appear  in  Section  IX,  and  conclusions  are 
given  in  Section  X. 


II.  THE  LAURITZEN- HOFFMAN  MODEL 

The  model  to  be  extended  is  one  version^’ ^  of  the  well-known  Lauritzen- 
Hoffman  (LH)  model  of  polymer  crystallization.  Our  description  of  this 
version  is  as  follows.  The  model  polymer  molecules  are  assumed  to  be 
infinitely  long  and  crystallize  on  an  existing  crystalline  face  or  substrate 
which  is  assumed  to  be  infinitely  wide  (i.e.  the  fact  that  its  width  is  finite 
is  ignored) .  A  sequence  of  length  i  of  polymer  segments  of  width  a  and 
thickness  b  as  well  as  the  volume  associated  with  that  sequence- -which  is 
taken  to  be  a  parallelepiped  of  length  i ,  width  a,  and  thickness  b--is 
designated  as  a  stem.  Only  stems  of  length  1  can  crystallize  on  an  existing 
face  of  length  I,  but  the  length  i,  the  lamellar  thickness,  can  vary  from 
crystal  to  crystal.  Any  sequence  of  length  I  of  segments  of  a  model  molecule 
can  be  placed  first  on  a  given  face  and,  upon  placement,  is  designated  as  the 
first  stem.  The  free  energy  of  formation  of  the  first  stem  is 

A^  -  A^q  -  A^  -  0  or  A^  -  2ab<7^  +  2bia  -  ablAf 

where  Af  >  0  Is  the  free  energy  of  fusion  per  unit  volume  at  a  temperature  T 
below  the  melting  point  of  the  model  polymer  (i.e.  of  a  crystal  of  very 
large  I)  and  Af  -  0  at  T  -  T’;  where  a  is  the  lateral  surface  free  energy  per 
unit  area  (i.e.  that  associated  with  the  surfaces  of  area  bl  and  a£  of  a 
stem) ;  and  where  o'  is  the  surface  free  energy  per  unit  area  associated  with 
the  cilium  that  protrudes  through  each  of  the  surfaces  of  area  ab  of  the  first 
stem.  Recently,^" ^  o'  has  been  assumed  to  be  zero;  generally,  one  can 
assume  that  0  <  o^  <  o&.  All  surface  free  energies  per  unit  area  in  the 
model  are  assumed  to  be  independent  of  T  and  I.  (See  Figure  2(a)  of  Reference 
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4  or  Figure  22  of  Reference  5.)  The  placement  of  each  subsequent  stem 
involves : 


1.  the  destruction  of  the  cilium  associated  with  one  of  the  surfaces  of 
area  ab  of  an  adjacent  stem  already  crystallized, 

2.  an  adjacent  reentry  and  the  formation  of  a  tight  fold  associated  with 
two  surfaces  of  area  ab,  and 

3.  the  formation  of  a  cilium  associated  with  the  remaining  surface  of 
area  ab  of  the  stem  being  placed. 

Only  adjacent  reentry  and  hence  only  tight  folding  is  incorporated  in  this 
version  of  the  model. 

The  free  energy  of  formation  of  the  i/th  stem  (i/  >  1)  is  therefore 


or 


A <f>v  -  ”  -abo^  +  2abag  +  abc^  -  abiAf 

A <f>  -  Atf  ,  -  2ab<x  -  abiAf  «  -E 
v  v-i.  e 


where  L4>v  is  the  free  energy  of  formation  of  a  group  of  v  stems  (relative  to 
A^q  ■  0)  and  where  a e  is  the  surface  free  energy  per  unit  area  associated  with 
half  of  a  fold.  Iteration  of  A^  -  A^ ^  •  E  (v  >  1)  gives 


A^  -  Atfx  -  (v-l)E 

-  2bia  +  2abOg  -  2aboe  +  »/ab(2ae-iAf) . 

In  order  that  stem  additions  subsequent  to  the  placement  of  the  first  stem  be 

thermodynamically  favorable,  i.e.  in  order  that  they  would  in  fact  occur,  one 

2a 

must  impose  the  constraint  -E  <  0  and  consequently  i  >  ^ - .  Stems  of  smaller 
length  are  unstable  and  disappear.  By  contrast,  A can  be  positive,  zero. 
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or  negative;  E  >  0  guarantees  that  A 4  <  0  will  occur  for  finite  v.  Note  the 


sign  conventions  for  A^  and  E. 


III.  THE  KINETIC  TREATMENT  OF  THE  LAURITZEN-HOFFMAN  MODEL 

Our  description  of  the  kinetic  treatment^’ ^  of  the  LH  model  is  as 
follows.  The  following  assumptions  are  made: 

1.  Assume  that  transition  state  theory  can  be  utilized  to  describe  the 
kinetics  of  the  LH  model  of  polymer  crystallization. 

2.  Assume  that  the  formation  (crystallization)  of  a  single  stem  is  an 
elementary  process  or  elementary  reaction;  that  the  destruction 
(melting)  of  a  single  stem  is  an  elementary  process  or  elementary 
reaction;  and  that  transition  state  theory  can  be  applied  to  these 
two  elementary  processes  with  a  single  transition  state  corresponding 
to  a  relative  free  energy  maximum  or  barrier  thus  occurring  between 
each  two  integral  values  of  v  on  a  plot  of  &4>v  vs.  v. 

3.  Assume,  that  only  one  stem  at  a  time  can  be  formed  or  destroyed. 

The  kinetic  problem  is  to  derive  an  expression  for  the  net  rate  S^i.T) 

at  which  stems  of  length  i.  (and  width  a)  pass  over  or  surmount  the  i/th  free 
energy  barrier  at  temperature  T.  The  problem  requires  consideration  of  the 
following  set  of  connected  elementary  reactions 

t°  h  h  h 

°  1  I  2  g  3  I  4  '  *  ' 

where  A  is  the  rate  constant  for  the  forward  reaction  v  -*  v+\  (v  >  1)  and  B  is 
that  for  the  reverse  reaction  i/+l  -*•  v  (v&l) ,  and  where  Aq  and  B-^  are  the 

-f  8 

analogous  rate  constants  for  the  u- 0  v-1  reactions.  Solution  of  this 

problem  in  the  steady- state  approximation  gives 


S„d.T) 


N0A0(A-B) 

A-B+B, 


S(i.T) 
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for  all  u,  whore  Nq  is  the  number  of  sites  or  locations  available  for  the 
placement  of  a  first  stem.  The  total  net  rate  at  which  stems  (i.e.  the  net 
rate  including  stems  of  all  possible  values  of  i)  pass  over  the  v th  barrier  at 
temperature  T  is  given,  for  all  v ,  by 


STotal(T)  “ 


S(i,T) 


i=J 


where  i^  is  the  smallest  allowed  value  of  i  which  satisfies  the  constraint 


i  > 


2_^e 
Af ' 


Note  that  i  is  a  discrete  variable- -  the  smallest  increment  in  i  that 


can  be  made  is  the  monomer  repeat  length  iu.  To  *ind  first  write  i  -  miu, 

2  a  2a /Af 

where  m  is  a  positive  integer.  Then  i  >  implies  m  >  —j -  ,  that  is,  m  is 

U2a  /Af 

greater  than  or  equal  to  the  smallest  integer  greater  than  — j -  ,  and 

"*u 

2a  /Af 

therefore,  i,  -  [l-*-INT(X)]i  ,  where  X  ■  — j -  and  INT(X)  designates  the 

i  -  u  j:u 

2oe 

integer  part  of  X.  Substituting  iu  ■*  into  the  expression  for  i^  gives  i^ 

-  [MEOll  fe]  .  To  .  good  approximation,  Ulffim  .  x  (i.e.  x  is 

2ae 

sufficiently  greater  than  1)  so  that  i ^  . 


Finally,  one  assumes  that 


S(i,T)  « 


i=3 


S(i,T)di;  and  the 


kinetically- determined  average  lamellar  thickness  is  then  given  by 

J®  i  S(i,T)  di 

i(T)  -  7^ - . 

J®  S(I,T)  di 
*1 


IV.  DETERMINATION  OF  THE  RATE  CONSTANTS 

To  obtain  expressions  for  Aq,  B^,  A  and  B,  one  must  first  determine 
expressions  for  the  free  energy  barriers  for  the  relevant  reactions  v  ^  v+1  (*/ 
>  0) .  Let  E^  be  the  free  energy  barrier  to  the  destruction  of  the  first  stem; 
then  A^j  +  E^  is  the  barrier  to  the  formation  of  the  first  stem  in  order  that 
(A^j+E^)  -  E-^  -  A^.  Let  E2  be  the  free  energy  barrier  to  the  formation  of 
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each  subsequent  stem;  then  E  +  E2  is  the  barrier  to  the  destruction  of  each 
such  stem  in  order  that  (E+E2)  -  E2  -  E.  Now,  one  does  not  know  the  free 
energy  barrier  to  the  formation  of  the  first  stem.  At  least,  one  does  know 
that  it  depends  on  what  length  i'  of  a  fully  adsorbed  stem  of  length  i 
actually  crystallizes  before  the  barrier  is  surmounted.  If  2'  -  0,  then  none 
of  the  free  energy  of  crystallization  (i.e.  -abiAf)  is  released  before  the 
barrier  is  surmounted,  and  clearly,  A0^  +  E^  -  2aba^+2b£a  and  E^  -  abiAf.  In 
general  then,  for  0  <£'<£, 

A0^+E^  -  2aba^  +  2bia  -  abi'Af  and  -  ab(i-i')Af. 

Since  £'  is  unknown,  a  parameter  0  ■  — j  with  0  ■&  0  <  1,  is  introduced  in  order 
that  all  possible  so-called  apportionments  of  the  free  energy  of  fusion  abiAf 
between  the  rate  constants  for  the  formation  and  destruction  of  a  first  stem 
(i.e.  for  the  forward  and  reverse  reactions  0^1)  can  be  considered.  Thus, 
A0-^  +  E^  -  2ab o'e  +  2bio  -  0abiAf  and  E^  -  (1-0;  abiAf. 

Note  that  the  greater  the  amount  0abiAf  of  the  free  energy  of  fusion  which  is 
in  fact  "apportioned"  (i.e.  the  greater  the  value  of  0  or  £'),  the  smaller  the 

value  of  both  A0,  +  E,  and  E,  (for  a  given  i  and  T) .  A  very  similar 
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interpretation  of  0  has  been  discussed  recently. 

Similarly,  for  each  subsequent  stem,  let  i"  (0  <  i"  <  i)  be  the  length  of 
a  fully  adsorbed  stem  which  actually  crystallizes  before  the  barrier  to  the 
formation  of  the  stem  is  surmounted.  Then 

E2  -  2abae  -  abi"Af  and  E  +  E2  -  ab(i-i")Af. 

Define  the  apportionment  parameter  $  ■  with  0  <  $  <  1  so  that 

E2  -  2ab(7g  -  $abiAf  and  E  +  Ej  -  (l-$)  abiAf. 


Finally,  utilizing  transition  state  theory, 
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A0~*t  e 


.T  -(A*,+Ei+A$)AT  -(A<£1+E1  )/kT 

^  e  11  -  0  e  11 


-E,  AT  -E,AT  -(E+E9)AT 

B-^  —  0  e  ;  A  —  e  ;  B  -  &  e 


where  a£  is  the  contribution  to  each  barrier  as  a  result  of  retardations  in 
the  transport  of  a  polymer  chain  through  the  liquid  to  the  substrate  or  vice 

B  A  ®i 

versa.  Note  that  -r  does  not  depend  on  \p  and  that  does  not  depend  on  ip  as 

A  Aq 

required. 


V.  THE  EXTENSION  OF  THE  IAURITZEN- HOFFMAN  MODEL 

As  implied  throughout  the  above  discussion,  the  application  of  transition 
state  theory  to  the  elementary  processes  of  single  stem  formation  and 

Q 

destruction  presumes  that  there  is  a  single  relative  free  energy  maximum  or 

barrier  between  each  two  integral  values  of  v  on  a  plot  of  A 4>v  vs.  v. 

Consequently,  Ac^+E^,  E^,  E2,  and  E+E2  should  never  be  negative  Clearly,  Et- 

(l-\p)  abiAf  and  E+E2  -  (l-$) abiAf  are  never  negative;  however,  the  expressions 

given  above  for  A^^+E^  and  Ej  can  be  negative.  In  fact,  E2,  for  example,  is 

2  O' 

negative  for  all  i  such  that  — —  <  i  for  a  given  Af,  \p,  and  a  .  We  propose  to 

0Af 

extend  the  LH  model  by  incorporating  into  the  model  the  assumption  that  free 
energy  barriers  cannot  be  negative,  i.e.  only  apportionments  of  the  free 
energy  of  fusion  which  result  in  a  nonnegative  barrier  will  be  allowed. 

In  order  to  incorporate  this  constraint  into  the  model,  first  note  that 
A^+E^  -  2ab a'e  +  2bicr  -  V’abiAf  is  never  negative  when  A<£^  is  positive  since 
then,  abiAf  <  2ab a'Q  +  2bio  always  holds  and  V’abiAf  <  2abcr^  +  2bio  follows. 
However,  when  A^  is  negative,  the  expression  2aba'&  +  2bicr  -  ^abiAf  can  be 
negative.  The  requirement  that  A^  +  E^  >  0  hold  when  A^  is  negative  implies, 
that  one  Is  not  allowed  to  apportion  all  of  the  free  energy  of  fusion  abiAf 
when  A<^  is  negative.  If  the  amount  ^abiAf  of  the  free  energy  of  fusion  which 
is  apportioned  were  to  exceed  2ab<7g  +  2bio,  then  A<^+E^  would  be  negative. 
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The  maximum  amount  which  can  be  apportioned  is  indeed  2aba^  +  2b£a,  and 
therefore  one  has,  when  A^  <  0, 

A^1+E]_  -  £ (2abo^  +  2b 2a) 

where  £  is  an  apportionment  parameter  with  0  <  £  <  1.  Using  (A<^+E^)  -E^  - 
A<^  or  E^  -  (A^+E^)  -  A^  gives 

E-,  -  C(2abo'  +  Ibla)  -  (2abc'  +  2bio-abiAf)  -  ab£Af- (l-<f )  (2abo'+  2blo) . 

Observe  that  the  requirement  that  A^+E^  >  0  holds  when  A^  is  negative  is 
equivalent  to  the  phys:  rally  realistic  requirement  that  the  barrier  E-^  to  the 
destruction  of  the  first  stem  cannot  be  smaller  than  the  free  energy  increase 
(-A<^)  that  occurs  upon  its  destruction.  Note  that  abiAf- (2ab<7M-2bio)  - 
-A<£^.  Also,  this  physically  realistic  requirement  implies  that  an  adsorbed 
first  stem  cannot  completely  crystallize  before  the  barrier  to  the  formation 
of  that  stem  is  surmounted,  i.e.  that  the  upper  limit  on  f  is  less  than  SL 
when  Atf^  is  negative.  This  upper  limit  on  i'  is  determined  later.  For  A<£^  > 
0,  the  expressions  A^  +  E^  -  2aba^  +  2bi<7  -  ^abiAf  and  E^  -  (l-^)abiAf  still 
hold  with  0  <  i/>  <  1  and  0  <  £'  <  i . 

At  this  point,  a  simple  change  of  variable  is  introduced  for  convenience. 
Define  A  -  l-£  with  0  <  A  <  1. 

Now  observe  that  although  the  free  energy  of  fusion  is  ab£Af  when  A<£^  is 
positive  or  negative,  the  free  energy  of  fusion  which  can  be  apportioned  is 
ab£Af  when  A <f>^  is  positive  but  is  (2abo^  +  2bio)  when  A^  is  negative.  Also, 
the  free  energy  of  fusion  that  is  in  fact  apportioned  is  ^abiAf  when  A<^  is 
positive,  but  is  A(2aba^  +  2b Ha)  when  A^  is  negative.  Clearly  then,  the 
fraction  of  the  free  energy  of  fusion  which  can  be  apportioned  that  is 
apportioned  is  rj>  when  A<£^  is  positive,  but  is  A  when  A^  is  negative.  If  we 
always  choose  the  same  value  for  A  and  \j>,  then  over  the  whole  range  of  values 
for  A^,  the  fraction  of  the  free  energy  of  fusion  which  can  be  apportioned, 
has  the  same  value.  Let  7  denote  any  particular  value  which  is  ctosen  for 
both  \f>  and  A ,  where  0  <  7  <  1 . 
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Note  that  equal  values  of  $  and  A  do  not  imply  the  same  value  of  £' 

0  f 

(except  when  A^  -  0  as  will  become  evident);  as  usual  i>  »  — j,  but  an 
expression  for  A  in  terms  of  i'  or  vice  versa  remains  to  be  obtained.  In  our 
approach,  then,  i'  depends  at  least  on  the  sign  of  A^ ,  and  yet  we  utilize 
only  one  parameter,  7- -the  fraction  of  the  free  energy  of  fusion  which  can  be 
apportioned  that  is  apportioned- -which  is  a  constant  over  the  whole  range  of 
values  for  A^. 

In  summary,  the  barriers  in  terms  of  the  apportionment  parameter  7  are 

A^  +  Ex  -  (l-7)(2abcr;  +  2bio) 

^  for  a^>2  -  0 

-  abiAf  -  7(2ab<7g  +  2bia) 

A <f>-^  +  -  2ab o'e  +  2bio  -  yabiAf 

I  for  Ac >  0 


Ex  -  (l-7)abfAf 

where  we  now  observe  that  (l-7)(2abo*  +  2bio)  -  2abo'  +  2bio  -  7abiAf  when  A^-, 

6  6  i. 

-  0,  i.e.  A$-^  +  is  a  continuous  function  of  i  and  Af  at  the  points  (i,Af) 

for  which  AcjJ^  -  0.  Note  that  the  greater  the  value  of  the  apportionment 

parameter  7,  the  smaller  the  value  of  both  A^  +  and  E^. 

An  expression  for  i'  is  not  needed  in  order  to  evaluate  ^^ota^(T)  and 

i(T).  However,  an  expression  for  V  in  terms  of  A  and  vice  versa  will  be 

derived  in  order  to  see  how  i'  depends  on  other  quantities  in  our  model. 

Given  A^  +  E2  -  (l-A)(2abo^  +  2bic)  for  A^  <  0,  or.e  can  first  find  >/>  when 

A<f>^  <  0  holds  in  terms  of  A  by  equating  the  expressions 

(1-A)  (2ab<7*  +  2bio)  -  2abo'  +  2bia  -  0abiAf 
0  6 


whence 

(-2ab  dg  +  2bfo’j 

*  "  A[  abiAf  J 


or 
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+  "  a(iaI  +  Sf)  • 

Clearly,  equating  these  expressions  and  expressing  \p  when  A<£^  <  0  in  terms  of 
A  is  valid  since  decreasing  2abcr^  +  2bio  by  an  amount  0ab2Af  must  be 
equivalent  to  decreasing  2abo^  +  2b 2a  by  A(2ab a'^  +  2bZa)  .  Note  that  the 

expression  +  aAf)  a^ways  less  than  one  when  A^  is  negative.  To  see 

this,  simply  observe  that  A<^  <  0  implies  2abo^  +  2bZo  <  abiAf,  and  then 

i ' 

divide  both  sides  of  this  inequality  by  abiAf.  But  ip  ■  —  f°r  all  values  of 
A^  so  that 

_  xif—  +  -2^1 

1  aAfJ  • 

Note  that  since  A  cannot  exceed  one,  the  largest  possible  value  of  £' ,  i.e. 

the  upper  limit  on  i',  is  i  for  A^  <  O’,  as  mentioned  previously, 

this  upper  limit  is  indeed  less  that  i  for  A^  <  0. 

For  completeness,  one  can  also  find  A  when  A<6^  >  0  holds  in  terms  of  ip  by 
equating  the  expressions 


(1-A)(2ab<7g  +  2bio)  -  2aba^  +  2bicr  -  ^abiAf 


whence 


_2  £  ) 
aAf  J 


Clearly,  equating  these  expressions  and  expressing  A  when  A<£^  >  0  in  terms  of 
ip  is  valid  since  decreasing  2abo^  +  2bia  by  an  amount  $abiAf  must  be 

0  * 

equivalent  to  decreasing  2abo'  +  2bla  by  A(2abo'  +  2bicr)  .  Here  again,  ip  -  — 


and  A  - 


1  f2ae  2c  ^ 

2^7 — -  .  Note  that  always  greater  than  one 

,  — —  "*■  — 
iAf  aAf 


when  A^  is  positive. 

In  summary,  then,  for  A^  <  0,  one  chooses  a  value  from  zero  to  one  for 

(2a'e  2  ol 

the  parameter  7,  whence  A  -  7,  and  then  calculates  ip  —  +  aAf  I  '  *’or  A^1 

>  0,  one  chooses  a  value  from  zero  to  one  for  the  parameter  7,  whence  0  -  7, 

V> 

and  then  calculates  A  -  - 2ctT'  '^iUS> 

(iAf  +  aAf] 
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A  -  7 


0  - 


^  for  <  0 


-  7 


A 


^  for  A^  >  0 


And,  for  all  A^ ,  one  can  calculate  i'  from  i'  -  ypt  or  from  2'-Ai^-|  +  j  • 
Incidentally,  the  constraint  2abOg  +  2bio  -  i^abiAf  >  0  combined  with  0  < 
ip  <  1  implies  that  the  inequality 


0  <  ip  <  the  smaller  of  1  and 


(ti  +  22_1 

(iAf  +  aAf  J 

must  be  satisfied,  and  clearly  our  theory  has  satisfied  it. 


Similarly,  the  constraint  abiAf  -  A(2abo^  +  2bicr)  >  0  combined  with  0  <  A  <  1 
implies  that  the  Inequality 

0  :£  A  <  the  smaller  of  1  and 


_L 


+  2S— 
[iAf  +  aAf 


must  be  satisfied,  and  clearly  our  theory  has  satisfied  It. 

The  approach  developed  above  can  readily  be  applied  to  incorporate  into 
the  model  the  constraint  that  E2  be  nonnegative.  Here,  E2  -  2aboe  -  $abiAf 
can  be  negative  when  E  Is  positive,  and  E  is  always  positive  (except  when  i  - 
2oe/Af,  which  gives  E  -  0) .  The  requirement  E2  >  0  implies  that  one  is  not 
allowed  to  apportion  all  of  the  free  energy  of  fusion  abiAf.  If  the  amount 
$abi Af  which  is  apportioned  were  to  exceed  2aboe ,  then  E2  would  be  negative . 
Therefore,  one  has  E2  -  *?2abcre  where  q  is  an  apportionment  parameter  with  0  < 
r)  <  1.  And  E+E-,  -  -2ab <7  +  abiAf  +  »72abc7  -  abiAf  -  (l-^)2abo  For 

4  6  6  C 
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convenience,  make  the  change  of  variable  ^  »  l-»/  with  0  <  ^  <  1  so  that  for 
all  i  and  Af 


E2  "  (l-0)2ab^e  and  E  +  Ej  —  abiAf  -  02abcrg. 


Observe  that  the  barrier  E  +  E2  to  the  destruction  of  the  second  and  each 
subsequent  stem  cannot  be  smaller  than  the  free  energy  increase  E  that  occurs 
upon  its  destruction,  which  implies  that  an  adsorbed  second  or  subsequent  stem 
cannot  completely  crystallize  before  the  barrier  to  the  formation  of  that  stem 
is  surmounted,  i.e.  that  the  upper  limit,  determined  below,  on  1"  is  less  than 
i. 

To  find  an  expression  for  £”  in  terms  of  8,  one  first  finds  $  in  terms  of 
8  by  equating  the  expressions  for  E2,  i.e. 

(l-0)2ab<r  -  2abo  -  #abiAf 

6  c 

whence 

iAf  • 


Clearly,  equating  these  expressions  and  expressing  $  in  terms  of  8  is  valid 
since  decreasing  2aboe  by  an  amount  $abiAf  must  be  equivalent  to  decreasing 
2abaft  by  02aboe.  Note  that  the  constraint  2ab<7g  -  $abiAf  >  0  implies  that  the 
A  2oe 

inequality  0  <  $  <  must  be  satisfied;  since  0  <  8  <  1  holds,  we  have 

2oe 

indeed  satisfied  this  inequality.  Also  note  that  -j£f  always  less  than  or 

2ae 

equal  to  one  since  £  >  has  been  established.  (Incidentally,  2abae  - 

a  2o  A»Af 

$abiAf  >  0  does  not  imply  constraints  £  <  ~ —  ,  Af  <  — E  ,  or  u  >  . ) 

yAf  e  c 


a  i"  e 

Finally,  recalling  that  V  ”  J~  and  substituting  above  gives  i"  - 

In  the  special  case  7  -  8  -  0,  our  mode^  reduces  to  the  case  if>  -  $  -  0  of 

the  LH  model  which  permits  negative  barriers  for  nonzero  rl>. 
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VI.  DETERMINATION  OF  THE  SIGN  OF  A^ 

At  this  point,  one  needs  to  determine  when  A^  is  positive,  zero,  and 
negative.  Now  A^  -  2aba^  +  2bia  -  abiAf  >  0  implies  b4(2a-aAf)  >  -2abcr^;  and 
there  are  three  cases  to  consider. 

Or,  -2abo' 

Case  (a):  2a  -  aAf  >  0  or  Af  <  Then  the  inequality  4  > 

is  always  satisfied  since  4  is  always  greater  than  zero,  and  hence  A^  >  0 

holds . 

Case  (b) :  2a  -  aAf  —  0  or  Af  —  Then  A^  —  2aba^,  which  is  always 

positive  or  zero  depending  on  o^. 

Thus,  combining  cases  (a)  and  (b) ,  we  have  A^  £  0  for  all  4  when  Af  < 

~,  where  a <f>,  -  0  when  both  a'  -  0  and  Af  - 
a  1  6  a 


Case  (c):  2 a  -  aAf  <  0  or  Af  > 


Then  A$^  >  0  implies  -b4(aAf-2a)  > 


M- 


-2aba '  or  4  <  — - 

e  -i.  2a_ 

1  aAf 


iQ.  Thus,  when  Af  >  A^  >  0  holds  for  4  <  iQ,  and 


A^  <  0  holds  for  4  2:  4q.  (Observe  that  as  Af  ^  from  values  greater  than 
^e:,  4Q  -  *>.)  There  is,  however,  one  further  condition  to  consider  here. 


2o  2a 

Recall  that  4  >  has  been  established.  If  4q  <  holds,  then  4  >  4q 

holds  and  consequently  A^  <  0  would  hold  for  all  4.  To  determine  when  4q  < 


2oL 

“Af  2cre  2  a 

s,  simply  write  - <  ~&£>  and  noting  that  <  1,  rearrange  this 

X"  aAf 


inequality  to  get  < 


—  "  .  Now,  if  og  5  o^,  this  inequality  would  be 


2ae 

<  0,  which  is  never  satisfied;  hence  4q  <  never  occurs  when  ag  <  a^.  If 
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2a* 

ae  > 

ae  *  i0  <  Af 

a 

>* 

f/v) .  * 

l<7e 

VII. 

EXPRESSIONS 

occurs  when  Af  > 


2a  f,  ge..J 
a  {oe-°'e  ] 


Thus,  if  and  ~  <  Af 

0  holds  for  £  >  £g,  but  for  Af 


(1) 

(2) 

(2) 


If  <r  <  o',  our  model  with 
e  e 

A^i+Ef  -  2abcr^+2b£a-7ab£Af 
A^+E^  -  2aboM-2b£a-7ab£Af 

A^^+E^  —  (I-7) (2abOg+2bIa) 


no  negative  barriers  has 

for  Af  < 

a 

for  Af  >  f2  and  £  < 
a 

for  Af  >  —  and  £  > 
a 


£ 


£ 


0 

0 


and  if  ag  >  a’e  , 

(1)  A<£}+E^  -  2aber^+2b£a-7ab£Af 

(2)  A^j+E-j^  -  2aba^+2b£o-7ab£Af 

(2)  A^1+E1  -  (l-7)(2abo;+2b£a) 

(3)  A^+Ej^  -  (l-7)(2aba^2b£a) 


for  Af  <  “■ 

cL 


and  £<£g 


and  £>£q 


The  purpose  of  categories  (1) ,  (2) ,  and  (3)  will  be  seen  shortly. 

When  A^+E-l  -  2abo'  +  2b£cr  -  7ab£Af,  Ex  -  (l-7)ab£Af,  which  we  call  Case 

I. 
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When  A<£1+E1  -  (l-7)(2ab<7^  +  2b la),  Ej_  -  abiAf  -  7<2ab<7^  +  2bio)  ,  which  we 
call  Case  II. 

One  always  has 


E2  -  (l-0)2ab<7£ 


E+E2  - 


■2aba  +  abiAf  +  E0  -  abiAf  -  82aba 

6  4  c 


Also, 


S(i.T) 


b 

'  !_  a  +  !i 

1  A  A 


B  -E/kT  B1  -(E1-E2)AT  j  b  -(A^+E^/kT 

where  t  —  e  '  ,  t~  -  e  ,  and  Aq  -  /3  e 


2aba ' 


Abbreviate  c' 


abAf 

kT 


ac 

i. 


2bo 

kT 


kT  * 

■  and  kT 


2abo£ 

~ kf“ 


2a 

aAf 


and  recall  i,  - 


£^e 

Af 


-  -c  +  q  i.  For  Case  I, 

*1 


Then 


A*i+Ei  .  2abg4  2bio  .  laMAf  .  c,  +  _c_ 

kT  kT  +  kT  kT  c  +  i, 


-7)i 


(1.7>1 .  <i.,)« 


For  Case  II, 


A«£1+E1  (I-7)  (2abo^+2big) 


<l-7)c'  +  7s  (l-7)i 


kT 


kT 


E, -E0  abiAf-7(2abc7'+2bio)  (l-0)2aba 

“kT"2 - V8 - ki^  -  if  a-«nr)<-y='-0-*)c 

For  Case  I, 

^Nne*c'  e-(«-7)c i/ix  (1.ec  e-c i/ij.) 

SI(t'T)  '  l-ec  ."=^1  -  e<l-*>c 

For  Case  II, 

fNfte-*1-*)®'  e-(1-7)aci/il  (l-ece-ci/£l) 

"  i-ece'ci/'el  +  e(1'*)c  e7C'  e’(1*a7 )ci/h 


Sn(i,T) 
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For  any  Af  in  category  (1),  then, 


slical(T>  '  ij-  f  dl  lnd  *' 


J.  1  S.(i,T)d  l 
l  _ 

SjCi.TJdi 


For  any  Af  in  category  (2) , 


stoLi<t>  -  i*r  f‘°  f.  SII<1^> 

u  JJCi  u  *0 

ft0  iST(i,T)di  +  H  i  S-r-r  (i,T)di 

j(2)m  .  Jl _ l _ tQ _ 11 

s^i.Ddi  +  Ji  Sjjd.Ud! 


For  any  Af  in  category  (3), 


^iSnd.Ddi 


m  i  r  n't  Jii  ‘'II'***/— 

STotal<T>  ‘ij,  SIl(I.T)di  and  i(3)(T)  -  -± - 

U  1  J i  Sjj >T)di 

For  purposes  of  comparison,  the  LH  model  which  permits  negative  barriers 
has,  for  all  t  and  Af, 

A^+E-^  -  2aba^  +  2bio  -  ^abiAf  and  E2  -  2aboe  -  $abiAf 


so  that 


E,  -E, 


-  (1-^+0)  i  -  c 


CLH'l  e'cV(a'^)ci/il  (l-ece'c'?/il) 

g\A-n^^  j)  —  - U - - - 

l-ec  e'^l  +  ec  e-d-’+4)c*/*i 

fS(LH)(i,T)di 

S&Sl*1-1)  -  i1  j!  S(LH)d,T)dl  and  l(Ut)(T)  -  1  - 

U  *1  J  .  S(LH)(^,T)di 

*1 

As  is  the  case  in  the  LH  model,  our  model  has  two  parameters.  The  most 


logical  choice  for  8  is  8  -  7;  however,  even  with  8  -  7,  our  integrals  cannot 
be  evaluated  analytically.  There  seems  to  be  no  special  case  (other  than  8  - 
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7-0)  for  which  they  could  be  evaluated  analytically.  At  this  point  then,  we 
proceed  without  setting  6  -  7. 


VIII.  EVALUATION  OF  THE  STotal(T)  AND  2(T)--THE  VARIABLE  TRANSFORMATIONS  FOR 
THE  NUMERICAL  INTEGRATIONS 

The  required  numerical  integrations  were  easily  performed  interactively 

10 

on  the  VAX  using  the  IMSL  subroutine  DQDAGS.  Integrals  to  be  evaluated  using 
DQDAGS  cannot  have  an  infinite  limit  of  integration.  One  way  to  proceed 
before  using  DQDAGS  is  to  make  a  change  of  integration  variable.  Although 
DQDAGS  can  integrate  functions  with  endpoints  singularities  (when  the 
endpoints  are  finite) ,  a  change  of  variable  which  results  in  a  transformed 
integrand  which  is  bounded  at  all  points  including  the  finite  e  dpoints  In  the 
new  range  of  integration,  is  preferable  to  a  change  of  variable  which  yields 
an  improper  integral  albeit  with  finite  integration  limits.  For  each  of  the 

integrals  appearing  in  sxotal^T)’  STotal^  '  and  STotal(T)  ’  a  variable 
transformation  which  resulted  in  a  proper  integral  was  in  fact  found.  The 
same  transformations  did  not  transform  the  corresponding  integrals  in  the 
numerators  of  2^(T),  2^(T),  and  2^(T)  into  proper  integrals;  however, 
the  transformed  integrands  were  of  the  form  (-In  x)f(x)  with  the  singularity 
resulting  only  from  the  factor  In  x  as  x  ■*  0.  This  endpoint  singularity  could 
be  handled  by  DQDAGS . 

Consider  first  the  integral  In  Sjotal^'  var*abl€  transformation 

consists  of  defining 

x  -  e^-^c  e’ (l*7)ci/f  1  . 

Note  that  x(i  -*  <*>)  -  0;  the  constant  ev  ,  i.e.  the  i  -  independent  factor, 

is  chosen  so  that  x(i  -  l^)  -  1.  Solving  for  i  in  terms 

Ji 

of  x  gives  l  -  2-jfl  -  provided  7*1-  Then  di  -  -  (^)dx. 

q--t  1 

Furthermore,  1  -  e-(a"*)c  x1*7,  e'ci/*l  -  e'c  x1"*  ,  and 
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-  e'd‘7)c  x  so  that 


Slital<T>  ‘  I 


u 


(l-7)c 


J°  ,  tH 


(i)  a* 


Simplifying  gives 


4ical<T> 


1-a-T  +  .I1'1)'  X 

a-1  1 


e  e  h  r1  x^Tfi-x1-^ 
“  i  (l-r)c  J0  _L. 

l-x1*7  +  e-<^-7)< 


dx 


This  is  one  of  the  integrals  that  was  evaluated  numerically  by  DQDAGS. 
Designate  the  integrand  above  as  f^(x) .  Using  the  same  variable 
transformation  to  evaluate  the  numerator  of  i^^(T)  gives 

Jo*!  ,  h  Jp(-ln  x)fx(x)dx 

"  "1  (l-7)c  rl 

Jo  fl< 


2(D(T) 


ll  fjCx) 


dx 


(x)  dx 


ro 

Next,  using  the  same  transformation  on  the  integral  Sj(i,T)  di 

appearing  in  4oLi<t>  elves 

,ta  e-c'e-(»-7)0f 

Sjd.TJdi  -  ^Nq  - k 


(l-l)c  Ix0fl<x>  * 


where 


Xq  -  x(i-£0)  -  e(1'7)c  e*(1'7)ci0/,ei  -  e(1'7)c  e-0-“7>cV(l-“) 


with  ip  ■  2c^/(l-a)Af  as  defined  previously. 

Similarly,  the  integral  iS][(i,T)di  appearing  in  £*2)(T)  becomes 


r*0  ^N0e'c'e'(a-7)c  £ 

JiJ  isi(i-T)  di  -  P . Ti-— 


1  {*i  Ix0fi<x>dx  +  rrfe  lx<-ln  x)ti<x)d>t} 


A  different  transformation  is  made  on  the  integral  J^Sjj(i,T)  di  also 
appearing  in  s£2£al(T).  Here,  define 

x  -  e(1*7)(c*c,)  B-(l'l)ac £/tl 
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Again  x(i  -♦<*>)  -  0;  the  constant  e^”7^cc  ^  i.r  chosen  so  that  x(i-iQ)  -  Xq  , 
which  is  given  above.  Solving  for  2  gives  2  -  ^  [c-c'-  ]  provided  7  * 

- —  (£)  dx.  Furthermore,  e- (l-7>aci/i1  „e'(l-7) (c-c' )  x 


1.  Then  di  -  -  (1.7)qc“" 

1 


e-*£/h  -  e-(c-c')/a  x(l-7)a  _  and  e* <1-07)0*/^  .  e 
Substituting  gives 

.-d-Tic'e-a-irK0-0')^ 


(c-c,)(l-Q7)  I.-.Q.7.. 


X 


d-7)a 


|i0sn(i,T)di“ 


r 

j  0 


1 


xl l-ece 


(l-7)ac 

HI 


,(l-7)«  1 _ ^ 

i-e°  7F?  j  x<rfe+  ea-*><=  .-re- 


^Nq  e'(1"r)ci1 

r  xo 

l-ece 

r  « *1 

c-c# 

1 

L  <*  J 

*d-7)« 

(1-7)qc 

0  l-ece  1 

c-c' 

l  l 

r«?-7)c  ec  e'| 

c-c' 

I  l-«7 _ 

l  « 

1  x(l-7)«  +  , 

L  <*  J 

1  xd-7)a 

-= -  (i)  dx 

1-cnr  vx' 

,(l-7)<* 


dx 


Designate  the  integrand  above  as  f2(x).  Similarly,  the  integral 
iSjjCi.TJdi  appearing  in  2^(T)  becomes 


r  ^N0e'(1*7)ci 

Jif  Il(I.T)di-  (1.7)ac 


1  *lj?  f2(x)dx+  (lW)ac  Io°('ln  x)f2(x)dx} 


Thus , 


s<2>  (T) 

bTotalu' 


(0 Nn  e'c'e'(a'7)ci,  rl  )  f/9N0  e-d-7)Cj»  (Xq  ^ 

-  (~T  (T^  Jx0fl(x)dX)  +  (  iu  (l-7)«c  Jo  f2(x)dxJ 


and 
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J(2)(T)  - 


iSI('e,T)di)  +  (/3N0  |ioiSll(i 


c(2) 

sTotal 


(T) 


/O') 

with  the  appropriate  expressions  for  the  integrals  and  sx<>tal^  to  de 
substituted  above. 

Finally,  consider  the  integral  in  The  variable  transformation 

to  be  made  on  this  integral  is 


x  m  e(l*7)QC  e-(l*7)ac^/^i 


Again  x(i  -»  «)  -  0  and  the  constant  e^*"^ac  is  chosen  so  that  x(i-i^)  -  1. 
Solving  for  i  gives  i  -  j^l-  ( 1 1  ^ac  ]  Provlded  7  »*  1.  Then  d £  - 

-  . (1)  dx.  Furthermore,  e‘ d-’7)*c*/*i  -  e-(l-7>«c  x  e*ci/il  - 

(l-7)ac  'x' 


_Lzsol 


,’c  x(l-7 )«  >  and  e-(l-a7)ci/i1  .  e-(l*a7)c  x(l-7>«  so  that 


s(3)  {T)  3  ! 

iTotal(-i;  i 


-(1-7) (c'+ac) » 

_ fj 

(l-7)ac 


ft 


(l-n(1~7)°) 


■i-Q,7, 


dx 


1.x(l-7)a  +  e-*c  e7(c'+txc)  x(l-7)a 


Designate  the  integrand  above  as  fj(x). 
evaluate  the  numerator  of  2(3)(T)  gives 

j(3)<t>  -  ^ + 


Using  the  same  transformation  to 


In  x)fj(x)dx 


^  fj(x)dx 


IX.  RESULTS  AND  DISCUSSION 

A  VAX  FORTRAN  program  was  written  to  evaluate  the  required  mathematical 
expressions.  All  calculations  were  done  double  precision  using  the  model 
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parameter  values  given  in  Figure  3  of  Reference  4;  namely,  a  ■  b  «  5  x  10 
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cm,  o  -  10  erg/cm^,  o&  -  100  erg/cm^ ,  -  500  K,  Ah  -  3  x  10^  ergs/cm^,  and 

o 

Af  -  (T^  -T)Ah/T^,  where  Ah  is  the  enthalpy  of  fusion  at  T  -  T  .  The  average 
lamellar  thickness  calculated  from  the  LH  model  is  independent  of  o^;  this  is 

true  for  our  model  only  for  Af  <  however.  Other  quantities  such  as 

11 

Sxotai(T)  do  depend  on  o'&  even  in  the  LH  model,  and  physically,  one  expects  0 

<  cr'  <  a  .  In  the  case  o'  -  0,  our  model  is  slightly  simpler,  for  then 
6  6  6 

A^  +  E^  -  2bio  -  yabiAf 
Zl  -  (l-y)abiAf 

A^  +  -  (l-y)2b 2o 

E^  —  abiAf  -  y2bicr 

Let  us  investigate  our  model  in  detail  for  the  case  a'e  -  0  first;  this  is  also 

1  2 

the  somewhat  arbitrary  choice  for  o^  made  for  the  calculations  ’  for  the  LH 

model.  For  the  values  of  a,  o,  Tm>  and  Ah  given  above,  the  temperature  T*  for 

which  Af  -  22  is  T*  -  433iK. 

a  3 

Given  the  parameter  values  above  and  now  with  the  choice  tf-7,  the 
calculated  average  lamellar  thickness  vs.  temperature  curves  (1  vs  T)  are 
plotted  in  Figure  1(a)  for  the  selected  values  of  7  -  0,  £,  and  (Results 

for  7  >  2  be  discussed  later.)  Some  of  the  data  used  to  construct  these 

plots  is  given  in  Table  I.  (For  Af  <  r2,  the  average  lamellar  thickness  is 
given  by  the  expression  for  2^(T)  given  previously  and  for  Af  >  by  the 
expression  for  2^(T)  also  given  previously.)  Clearly,  2  decreases 
monotonically  with  decreasing  T  in  agreement  with  typical  experimental 
behavior.  For  most  supercoolings,  the  magnitude  of  the  2  values  is  of  the 
order  of  25-125A,  which  is  quite  reasonable.  Note  that  at  least  for  all 
values  of  Af  >  2  at  a  given  T  increases  with  increasing  7.  Also,  the 

numerical  results  shown  in  Figure  1(a)  indicate  that  2  vs.  T  is  relatively 
insensitive  to  the  value  of  7. 


23. 


For  comparison,  we  have  reproduced  part  of  Figure  3(b)  of  Reference  1  as 
our  Figure  1(b),  which  shows  the  1H  model  2  vs,  T  curves  with  $  -  ip  for  the 
selected  values  of  ^  -  0,  £,  and  Some  of  the  data  which  we  calculated 

in  order  to  construct  these  plots  is  given  in  Table  II.  The  LH  model  i>  -  0 
curve  is  identical  to  our  7  -  0  curve.  For  Af  <  each  of  the  LH  model  "\p 
curves"  is  qualitatively  similar  but  not  quantitatively  identical  to  its 
corresponding  "7  curve"  presented  in  Figure  1(a).  Recall  that  the 
quantitative  difference  arises  from  the  fact  that  the  barrier  E2  has  been 
constrained  to  be  nonnegative,  i.e.  E2  ”  (l-0)2aboe.  For  Af  >  ,  however, 

the  LH  model  ip  curves  are  in  marked  contrast  to  the  7  curves;  in  particular, 
for  each  ip  curve,  2  approaches  infinity  asymptotically  as  Af  approaches  ^ . 
This  is  the  behavior  which  is  known  as  the  62  catastrophe. 

One  point  is  worth  emphasizing  here;  namely  the  relationship  between  7 

0 1 

and  ip.  In  both  our  model  and  the  LH  model,  iP  ■  but  this  ratio  in  the  LH 
model  is  a  constant,  whereas  in  our  model 


*'{ 


,  rife  +  _22_1 

7  (iAf  aAf  J 
7 


A^  <  0 
A4>1  >  0 


For  the  case  a'  -  0,  this  becomes 
e 


-{ 


y  -2°. 
7  aAf 


> 


Af  > 


Af  <  ~ 
a 


Now,  for  any  given  ip,  say  ip.  ,  I  in  the  LH  model  Is  Infinite  for  all  Af  >  ~r~ 

J  Vja 

a 

and  for  all  Af  >  there  is  no  finite  value  of  2  for  any  ip  >  ip j  . 

Equivalently,  a  value  of  ip  >  ip^  is  not  possible  for  a  chain-folded  system  for 

all  Af  >  -r^2,  that  is,  high  values  of  ip  do  not  lead  to  chain- folded  polymer 
V>  j  a 

crystals  at  high  enough  supercooling  according  to  the  LH  model.  Experiment, 
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however,  gives  chain- folded  crystals  at  high  supercooling  with  an  average 

lamellar  thickness  that  decreases  monotonically  with  decreasing  temperature. 

As  we  have  seen,  our  one-parameter  (i.e.  7)  model  with  o'  —  0  does  reproduce 

this  high  supercooling  behavior.  And  yet,  high  values  of  i.e.  of  the  ratio 
i' 

are  not  associated  with  our  high-supercooling  chain-folded  systems.  To 

see  this,  first  introduce  the  dimensionless  quantity  x,  where  0  <  x  <  1.  Then 

2a  2 

for  any  Af  -  — ,  -  7  -  7X.  Since  7  cannot  exceed  one,  ^  In  our  model 

cannot  exceed  x.  for  any  Af  >  — where  x.  is  any  given  value  of  x.  But  this 
J  x  j a  J  '  ■ 

is  exactly  what  was  found  for  if>  in  the  LH  model,  i.e.  that  a  value  of  ^ 

greater  than  or  equal  to  is  not  possible  for  any  Af  >  .  Thus,  for  Af  > 

j  *ja 

our  model,  through  the  imposition  of  the  constraint  that  barriers  be 

nonnegative,  places  exactly  the  same  upper  limit,  ,  on  our  \p  that  is 

predicted  for  in  the  LH  model.  However,  for  Af  >  — ,  our  model,  unlike  the 

<£l 

LH  model,  predicts  1  vs,  T  in  qualitative  agreement  with  experiment. 

Th..s,  the  selected  calculations  done  for  our  model  indicate  that,  for  the 
case  a'e  -  0,  our  model  does  pot  exhibit  an  infinite  average  lamellar 
thickness.  Most  importantly,  our  model  predicts  1  vs.  T  curves  which  are 
monotonically  decreasing  with  decreasing  T  in  agreement  with  experiment.  That 
is,  we  have  successfully  extended  the  LH  model  to  higher  supercooling. 

Also,  this  success,  coupled  with  the  numerical  results  shown  in  Figure  1(a), 

significantly  increases  our  confidence  in  using  7  -  0  as  a  first  approximation 

7 

for  mathematical  convenient?  in  practice.  Finally,  our  results  show  that  the 
Si  catastrophe  of  the  LH  theory  is  related  to  the  failure  to  exclude  negative 
barriers  and  moreover  that  the  LH  approach  to  polymer  crystallization  is  in 
itself  valid  for  high  supercooling- -given  that  negative  barriers  are 
forbidden.  Prior  to  this  work,  the  LH  approach  had  always  been  described  as 
one  which  is  invalid  at  high  supercooling. 


One  set  of  results  with  8  *  7  is  presented  in  Table  III.  Here  we  see 
that  for  7  “  2  an<*  ^  che  calculated  1(T)  differ  only  slightly  from  the 

case  with  7  -  ^  and  $  - 
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Next,  we  investigated  our  model  for  a'^  *  0.  (Recall  that  2  for  the  LH 
model  is  independent  of  o’e  and  that  our  model  is  independent  of  a '  for  Af  < 
^.)  Using  the  same  values  for  a,  b,  a,  a g,  T^,  and  Ah  as  above  and  again 
with  0  -  7,  2  vs.  T  curves  for  a'&  -  0,  60,  100,  and  150  erg/cm  --each  with 
7  -  ^--are  plotted  together  in  Figure  2.  Some  of  the  r*  0  data  used  to 
construct  these  plots  is  given  in  Table  IV  (and  the  -  0  data  has  been  seen 
previously  in  Table  I).  From  Figure  2,  we  see  that  2  decreases  monotonically 
with  decreasing  T  for  0  <  as  well  as  for  a'e  -  0  and  that  2  vs .  T  is 

relatively  insensitive  to  the  value  of  a ^  <>  a&.  Thus  our  conclusions  made 
immediately  above  for  the  case  a'e  -  0  are  valid  when  0  <  s  ag.  For  - 

O 

150  erg/cm  ,  there  is  a  relative  minimum  in  2  vs.  T  near  T  -  405  K,  and  the 

curve  passes  through  a  small  and  "diffuse"  relative  maximum  at  a  lower 

temperature.  Recall  that  one  expects  0  <  o'  <  a so  that  with  a  - 

6  6  6 

2  2 

100  erg/cm  ,  a'&  -  150  erg/cm  may  not  be  realistic  but  is  examined  in  order  to 
explore  the  model  predictions  as  a  function  of  a '. 

The  relationship  between  7  and  ^  with  <*'e  *  0  is  worth  emphasizing  at  this 


point.  In  doing  so,  one  difference  between  the  cases  a'^  -  0  and  cr^  1*  0  will 
be  found;  namely,  i>  can  exceed  for  some  Af  £  -7^  when  cr'&  *  0.  To 

i' 

reiterate,  in  both  our  model  and  the  LH  model,  V>  “  J~,  hut  this  ratio  in  the 


LH  model  is  a  constant,  whereas  in  our  model 

f2a'_  -v 


!,T)  -  - 

k. 


+  22-. } 
+  aAfJ 


A^X(2,T)  <  0 
A^1(2,T)  2:  0 


where  the  notation  \6(2,T)  and  A<^(i,T)  emphasizes  here  the  dependence  of  r(>  and 
A^  on  l  and  T.  (The  T  dependence,  of  course,  enters  through  Af.)  Recalling 
the  conditions  which  govern  the  sign  of  A^  then  gives,  when  >  a'e 


26. 


[I a_k  .  _2a] 

7  [iAf  +  aAfj 


*<i.T)  -  \ 


{ 


for  all  i  when  Af  > 


2a 


(^V) 

v  P  P-' 


for  i  >  £q  when 


2a 


<  Af  < 


< 


a 


fe) 


{ : 


la 


for  i  <  20  when  f2  <  Af  <  f2  [~r] 

v  a  py 


or  all  i  when  Af  5 
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and  when  a  <  a' 
e  e 


*(i.T) 


-1 


-  ffft  +  2s_] 

7  (.iAf  + 


aAfj 


for  2  >  2n  when  Af  >  ^ 
0  a 


{ 


for  2  s  iQ  when  Af  >  a 


2a 


for  all  2  when  Af  < 


2s 


where  2q 


Furthermore,  on  an  1  vs.  T  curve,  one  has 


A^d.T)  5  0 
A^d.T)  >  0 


where  the  conditions  which  govern  the  sign  of  A^(2,T)  are  those  given  above 
for  A^i.T)  but  with  £  replaced  by  2.  Therefore,  the  temperature  TQ  of  a 
point  (2q,  Tq)  on  an  2  vs.  T  curve  and  at  which  A^(2,T)  -  A^1(2q,Tq)  -  0  is 
the  solution  to  the  following  non-linear  algebraic  equation  in  the  one  unknown 
T: 

2(2)(T)  -  £q 
or 
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Ji1  iSIdi+Ji0iSIIdi 
SIdi+  ^0SIIdi 


2a' 


_4L 


aAf 


» «e  > 


T0  will  correspond  to  a  value  of  Af  in  the  range  jj2  <  Af  < 


%o_ 

a 


Af  > 


e  eJ 

Zs. 

a‘ 


but  if  a  <  a' 
e 


,  2q  will  correspond  to  a  value  of  Af  in  the  range 


Rather  than  attempt  to  solve  the  above  equation  iteratively,  one  simply 
plots  the  left-hand  side  2^(T)  vs.  T  and  the  right-hand  side  ig(T)  vs.  T  on 
the  same  graph,  and  Tq  is  given  by  a  point  of  intersection  of  the  two  curves. 


Note  that  as  Af  approaches  —  from  values  greater  than  ^2, 


Iq  approaches 


2a 


infinity  and  that  ig  decreases  monotonically  with  decreasing  T  for  Af  >  &  . 

For  each  of  the  2  vs.  T  curves  with  r*  0  ,  we  found  one  point  of 
intersection  (2q,  Tq) ,  which  is  designated  on  each  curve  by  an  open  circle. 

We  also  found  that  2^(T)  >  holds  when  T  <  Tq  and  that  2^(T)  <  2q  holds 
when  T  >  Tq.  Thus,  A^(2,T)  <  0  holds  for  T  <  Tg  and  A^(2,T)  >  0  holds  for  T  > 


lo- 


Our  final  result  is  that,  on  an  2  vs.  T  curve, 

r2a' 


r  »  P + 

-  4  [2Af  aAfJ 


0  <  T  <  T, 


To  *  T  <  t; 


Note  that  if  the  dimensionless  quantity  x,  0  <  x  <  1,  is  again  introduced  by 

/•aal 

writing  Af  -  fr,  then  ^ 


so  that,  unlike  the  case  a'&  -  0,  can 


exceed  Xj  for  some  Af  >  ^ -,  where  x^  is  any  given  value  of  x. 

Now,  upon  proceeding  to  consider  results  for  7  >  ^  ,  our  basic 
conclusions- -especially  the  fact  that  we  have  removed  the  52  catastrophe  at 
high  supercooling- -remain  intact;  however,  we  do  not  obtain  2  vs.  T  curves 
which  are  monotonically  decreasing  for  all  T  when  7  is  "sufficiently"  large. 
Using  the  same  values  for  a,  b,  a,  ag,  Tn° ,  and  Ah  as  previously  and  again 
with  5-7  and  o'&  -  0,  the  calculated  2  vs.  T  curves  for  the  selected  values 
of  7  -  .90,  and  .95  are  plotted  in  Figure  3(a),  and  the  curve  for  7  - 
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.99  appears  in  Figure  4.  Some  of  the  data  used  to  construct  these  plots  is 
given  in  Table  V.  The  effect  of  7  on  2  as  a  function  of  T  is  readily 
apparent.  First,  the  curve  for  7  -  2  aPPears  on  closer  examination,  to 
exhibit  a  discontinuity  or  break  in  its  slope  at  the  temperature  T*  -  433^K 
for  which  Af  -  (This  statement  will  be  qualified  later.)  As  for  7  - 

O 

2  for  7  -  £,.9,  .95,  and  .99  does  decrease  with  decreasing  T  for  all  T  for 
which  Af  >  and  there  appears  to  be  a  break  in  the  slope  of  2  vs.  T  at  T  - 
T*.  Unlike  for  1  “  the  higher  7  curves  pass  through  a  relative  minimum  at 
a  temperature  for  which  Af  <  the  temperature  T*  at  which  this  minimum 
occurs  increases  with  7  (for  7  -  £,  it  occurs  between  T  -  440  and  433-jK  and  so 
can  hardly  be  seen  on  the  plot.)  Also,  over  the  interval  T  <  T*,  2  vs.  T  is 
at  a  relative  maximum  at  T  -  T*.  Finally,  note  that  2  vs.  T  curves  for  .99  < 

7  <  1  are  qualitatively  similar  to  the  7  -  .99  curve  and  do  not  exhibit  an 
infinite  average  lamellar  thickness.  The  numerical  integrations  in  the 
expressions  for  2^^(T)  and  2^^(T)  could  not  be  done  for  7  ■  1  as  a  result  of 
the  factor  (I-7)  appearing  in  various  denominators. 


For  comparison,  we  have  reproduced  part  of  Figure  3(b)  of  Reference  1  as 

our  Figure  3(b),  which  shows  the  LH  model  2  vs.  T  curves  with  $  -  i/>  for  the 

1  3 

selected  values  of  rf>  -  4.  .90,  and  .95.  Some  of  the  data  which  we 

calculated  in  order  to  construct  these  plots  is  given  in  Table  VI .  These  IH 
model  if>  curves  exhibit  the  Si  catastrophe  as  Af  approaches  J§,  as  do  all  LH 
curves  for  . 95  <  3  1.  The  curves  for  . 95  <  0  <  1  are  similar  to  the  xf>  -  .95 

curve;  since  integrations  can  be  done  analytically  in  the  LH  model  when  $  - 
2  vs.  T  for  i>  ■  1  was  able  to  be  obtained.^" 

Thus,  for  high  enough  7,  our  <7^-0  model  2  vs.  T  curves  appear  to  have  a 
break  in  slope  at  T  -  T*.  We  suspect  that  there  is  indeed  a  break  in  slope  at 


7 
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implies  that  ^  is  discontinuous  at  Af  -  however,  we  have  not  evaluated  ~~ 

at  Af  -  The  break  in  slope  is  apparently  indiscernible  up  to  7  values  of 

about  where  the  slope  of  2  vs.  T  has  the  same  sign  (positive)  regardless  of 

whether  the  point  Af  -  “  is  approached  from  values  of  Af  higher  or  lower  than 
r!<- 

As  7  increases,  however,  the  break  becomes  pronounced  with  the 

a 

concomitant  appearance  of  a  relative  maximum  in  1  at  T  -  I*  and  a  relative 

minimum  in  2  at  T  -  T*;  necessarily  then,  the  slope  of  2  vs.  T  as  Af 

approaches  ^2.  from  values  less  than  ^  becomes  negative.  We  will  refer  to 

this  undesirable  behavior,  manifest  at  high  values  of  7,  as  the  2  anomaly. 

Unlike  the  Si  catastrophe  in  the  LH  model,  the  relative  maximum  in  2  vs.  T,  as 

noted  above,  always  appears  at  Af  -  ^  for  all  values  of  7  given  that  o'e  -  0. 

Next,  we  consider  o'  *  0  for  high  values  of  7,  The  2  vs.  T  curves  for  o' 

2  3 

-  0,  60,  100,  and  150  ergs/cm  --each  with  7  "  4  --are  presented  in  Figure  5. 
The  curves  pass  through  a  common  relative  minimum  between  T  -  440  and  433^  K 
(for  which  Af  <  ) ,  and  then  each  curve  rises  and  passes  through  a  relative 

ft 

maximum,  that  maximum  being  relatively  higher  and  occurring  at  higher  Af  the 
larger  the  value  of  a • .  At  each  maximum,  there  would  appear  to  be  a  break  in 
the  slope  of  2  vs.  T.  Having  passed  through  its  maximum,  each  curve  decreases 
monotonically  with  decreasing  T  thereafter. 

One  should  be  careful  to  note  that  what  appears  to  be  a  break  in  the 

d2 ( (T\ 

slope  of  2  vs.  T  when  o'&  *  0  is  probably  not  a  break  in  slope;  - — L  should 

be  continuous  for  all  relevant  T.  Whether  a  break  in  the  slope  of  2  vs.  T 
occurs  at  Af  -  “  when  o'  *  0  as  was  presumed  true  for  o'  -  0  cannot  be 
determined  conclusively  from  the  appearance  of  the  graphs,  although  the  break 
appears  to  be  absent. 

Qualitatively  similar  2  vs.  T  curves  are  obtained  for  7-0.9  and  o^  -  0, 

0 

60,  100  and  150  ergs/cm  as  is  shown  in  Figure  6.  See  also  Table  VII.  The 
relative  maxima  are  higher  and  "sharper"  than  the  corresponding  7  “  4  curves, 
and  they  have  moved  to  higher  temperature.  For  7  -  0.99,  the  analogous 
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curves,  shown  in  Figure  7,  exhibit  1  values  which  are  unrealistically  large  as 
well  as  maxima  which  are  extremely  "sharp". 

Thus ,  from  the  graphs ,  we  see  that  the  1  anomaly  becomes  more  pronounced 
but  moves  to  higher  temperature  as  7  increases  for  a  fixed  nonzero  value  of 
o'  That  is,  although  the  relative  maximum  in  I  vs.  T  can  appear  at  some  Af  > 

A 

~  when  0^  is  nonzero,  the  maximum  becomes  less  pronounced  as  it  moves  to 
lower  temperature  upon  a  decrease  in  7.  Our  model,  then,  does  not  fail  at 
high  supercooling,  but  does  exhibit  anomalous  behavior  for  temperatures 
corresponding  to  values  of  Af  "just"  greater  and  "just"  less  than  ^ .  This 
undesirable  behavior  is  pronounced  for  large  values  of  7  and  is  more 
pronounced  for  larger  values  of  o'&  for  a  given  7. 

We  can  easily  rationalize  mathematically  how  our  calculated  2  vs.  T 
curves  can  rise  with  decreasing  T  for  some  Af  >  ^  when  o'  is  nonzero.  Recall 
that  the  expression  for  J^^(T),  namely 


contains  two  different  integrands  Sj(i,T)  and  Sjj(i.T).  Depending  on  o'e,  7, 
and  T,  the  contribution  of  the  integrals  involving  Sj(i,T)  to  1^(T)  may 
outweigh  the  contribution  of  the  integrals  involving  Sjj(i,T),  and  in  some 
cases,  our  calculations  show  that  to  a  very  good  approximation 


If0  iSI(i,T)d £ 

2^(T)  «  - - -  with  2q  approaching  infinity. 

Sj(i  ,T)df 

But  this  is  our  expression  for  2^(T)  for  the  interval  Af  <  and  the 
results  of  our  calculations  using  l^CT)  have  been  found  to  differ  little 


from  results  using  i^^T),  i.e.  the  LH  theory.  Not  unexpectedly  then, 
j(2)^)  can  increase  with  decreasing  T  for  some  Af  >  ^.  We  note  that  the 
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numerator  of  S^(i,T),  like  the  numerator  of  contains  the  factor  Aq 
-  e*c  e*hi(2a*7aAf)/kT^  t^e  form  Qf  which  has  been  associated  vith^increascs 
in  2  with  decreasing  T. 


X.  CONCLUSIONS 

We  have  constructed  a  model  of  polymer  crystallization  which  extends  the 
LH  theory  by  excluding  negative  free  energy  barriers ,  and  we  have  shown  that 
the  61  catastrophe  of  the  LH  theory  is  related  to  the  failure  to  exclude  these 
negative  barriers.  Our  results  show  that  the  new  model  is  more  consistent 
with  experimental  behavior  at  very  high  supercooling. 

Our  results  with  -  0  clearly  indicate  that  the  2  anomaly  in  our 
model* -and  in  part  the  61  catastrophe  of  the  LH  theory- -are  associated  with 
the  interval  Af  S  ^  and  are  thus  connected  to  the  expression  A^  +  -  2ab<7^ 

+  2b la  -  yabiAf.  The  2  anomaly  also  appears  to  be  connected  to  this 
expression  even  when  a'&  y*  0,  i.e.  even  when  the  maximum  in  2  vs.  T  occurs  at 
a  temperature  for  which  Af  exceeds  Although  high  values  of  7  and  \/>  are 

considered  unrealistic  as  has  been  elucidated^  recently,  however,  there  is  no 
guarantee  that  the  LH  theory  as  well  as  our  extension  of  it  has  not  failed  to 
incorporate  an  as  yet  unknown  constraint  or  feature  which  would  improve  the 
model  results  at  high  7  values.  For  example,  high  7  values  may  be 
unrealistic,  but  the  2  values  for  high  7  from  an  improved  model  may  simply  be 
unrealistically  large  but  nevertheless  monotonically  decreasing  with 
decreasing  T  for  all  T. 

In  conclusion,  we  hope  to  extend  our  modification  of  the  LH  approach  to 
polymer  crystallization  to  treat  the  interesting  systems  which  interact  with 
an  applied  electric  field. 
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FIGURE  CAPTIONS 


Figure  1(a).  Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 

for  7-O,  and  each  with  a’&  -  0  and  6  -  7.  See  Section 

IX  for  a,  b,  <7,  cx  ,  T®,  and  Ah  which  are  the  same  for  all  of 

6  IB 

the  figures.  At  T  -  433^  K  (i.e.  Af  -  ^)  ,  A^  -  0. 

For  T  >  433^  K,  A^  £  0  and  16-7  and  A  -  7^|J  • 

For  T  <  4332  ^1  ~  ®  and  ^  “  7(aAf )  anc*  A  —  7. 


Figure  1(b).  Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 
for  rf>  -  0,  £,  j,  and  xjr,  each  with  0  -  rf>,  reproduced  from  the 

Lauritzen-Hoffman  Model  (Reference  1) ;  plots  are  independent 

of  o'. 
e 


Figure  2.  Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 

for  o^  -  0,  60,  100,  and  150  ergs/cm^,  each  with  8  -  7  -  j. 
Each  open  cirle  designates  the  point  (Iq,  Tq)  at  which 
A^2(2,  T)  -  0.  For  T  2  Tq,  A^  st  0,  ^  -  7,  and 

.  y(ti  +  _2d 

W  “fl' 


..ablA.f _ 

2abtx'  +  2bl<x 
e 


-) 


For  Ti  Tn,  A^1  <  0,  \j> 


and  A  -  7. 


Figure  3(a).  Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 
for  7  -  0.90,  and  0.95,  each  with  o'&  -  0  and  8  -  7. 

As  in  Figure  1(a),  A^  -  0  at  T  -  433^  K. 

Figure  3(b).  Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 

for  V>  "  2'  4'  and  .95,  each  with  0  -  reproduced  from 

the  Lauritzen-Hoffman  Model  (Reference  1) ;  plots  are 
independent  of  0 

Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 
for  9  -  7  -  .99  and  o ^  -  0.  As  in  Figure  1(a), 

A^1  -  0  at  T  -  433^  K. 


Figure  4. 
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Figure  5. 


Figure  6. 


Figure  7 . 


Plots  of  Average  Lamellar  Thickness  (A)  vs .  Temperature  (K) 
for  cr^-O,  60,  100,  and  150  ergs/cm2,  each  with  8  -  7  -  £. 

As  in  Figure  2,  each  open  circle  identifies  the  temperature  Tq 
at  which  A^^(5,T)  -  0. 

Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 
for  <7^-  0,  60,  100,  and  150  ergs/cm2,  each  with  8  -  7  -  .90. 
As  in  Figure  2,  each  open  circle  identifies  the  temperature  Tq 
at  which  A^(l,  T)  -  0. 

Plots  of  Average  Lamellar  Thickness  (A)  vs.  Temperature  (K) 
for  o'  -  0,  60,  100  and  150  ergs/cm2,  each  with  8  -  7  -  .99, 
For  o'B  -  0,  60,  100,  and  150  ergs/cm2,  TQ  -  433^  K,  432.2  K, 
432.1  K,  and  432.0  K,  respectively.  As  in  Figure  2,  Tq  is  the 
temperature  at  which  A^(1,T)  -  0. 


TABLE  CAPTIONS 


Table  I.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 
for  7-O  and  for  7  -  each  with  0^  -  0  and  0  -  7.  See  Figure 

O 

1(a).  See  Section  IX  for  a,  b,  a,  oe,  Tm,  and  Ah,  which  are  the 
same  for  all  of  the  tables. 

Table  II.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 
for  ij>  -  2  with  0  -  \f>,  reproduced  from  the  Lauritzen-Hoffman 
(LH)  Model  (Reference  4);  data  is  independent  of  o’.  See  Figure 
1(b). 

Table  III,  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 
for  7  “  2*  0  “  1,  and  o'e  -  0. 

Table  IV.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 

2  1 

for  o'e  -  60,  100,  and  150  ergs/cm  ,  each  with  0  -  7  -  See 
Figure  2. 

Table  V.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 
for  7  -  .90  with  0^  -  0  an<|  0-7.  See  Figure  3(a). 

Table  VI.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 

for  ^  -  .90  with  $  -  reproduced  from  the  Lauritzen-Hoffman  (LH) 
Model  (Reference  4);  data  is  independent  of  0^.  See  Figure  3(b). 

Table  VII.  Average  Lamellar  Thickness  (A)  as  a  function  of  Temperature  (K) 
for  o'e  -  60,  100,  and  150  ergs/cm2,  each  with  0  -  7  -  .90. 

See  Figure  6. 
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Table  I.  Average  Lamellar  Thickness  l R) 
vs.  Temperature  (K) 

TEMP.  (K)  Psi«Gamma«0  Gamma- 1/2 


485.000 

480.000 

475.000 

470.000 

465.000 

460.000 

455.000 

450.000 

445.000 

440.000 

435.000 

430.000 

425.000 

420.000 

415.000 

410.000 

405.000 

400.000 

395.000 

390.000 

385.000 

380.000 

375.000 

370.000 

365.000 

360.000 

355.000 

350.000 

345.000 

340.000 

335.000 

330.000 

325.000 

320.000 

315.000 

310.000 

305.000 

300.000 

295.000 

290.000 

285.000 

280.000 

275.000 

270.000 

265.000 

260.000 

255.000 

250.000 

245.000 

240.000 

235.000 


234.383 

178.390 

144.660 

122.074 

105.867 

93.652 

84.105 
76.429 
70.115 
64.826 
60.328 
56.451 
53.072 
50.100 
47.463 

45.105 
42.984 
41.064 
39.316 
37.718 
36.251 
34.897 
33.644 
32.480 
31.396 
30.382 
29.433 
28.540 
27.700 
26.907 
26.157 
25.446 
24.772 
24.130 
23.518 
22.935 
22.377 
21.843 
21.332 
20.841 
20.369 
19.915 
19.479 
19.057 
18.651 

18.258 
1  7.878 
17.511 
17.155 
16.809 
16.475 


235.303 

179.781 

146.556 

124.507 

108.867 

97.253 

88.342 

81.344 

75.762 

71.267 

67.641 

63.528 

59.481 

55.988 

52.941 

50.259 

47.877 

45.744 

43.821 

42.077 

40.484 
39.023 
37.676 
36.429 
35.270 
34.188 
33.176 

32.225 
31.329 

30.484 
29.683 
28.924 
28.201 
27.513 
26.855 

26.226 
25.624 
25.045 
24.489 
23.953 
23.437 
22.938 
22.456 
21.990 
21.537 
21.099 
20.673 

20.259 
19.856 
19.463 
19.081 


Table  II.  Average  Lamellar  Thickness  (£)  vs.  Temperature  (K) 


TEMP  (K)  LH  Psi“1/2 


485.000 

235.785 

480.000 

180.224 

475.000 

146.926 

470.000 

124.780 

465.000 

109.027 

460.000 

97.290 

455.000 

88.251 

450.000 

81.124 

445.000 

75.412 

440.000 

70.789 

435.000 

67.037 

430.000 

64.009 

425.000 

61.610 

420.000 

59.786 

415.000 

58.519 

410.000 

57.832 

405.000 

57.800 

400.000 

58.577 

395.000 

60.458 

390.000 

64.019 

385.000 

70.494 

380.000 

82.999 

375-000 

112.171 

370.000 

232.547 

365.000 

oo 

Table  III.  Average  Lamellar  Thickness  (%.)  vs.  Temperature  (K) 


TEMP.  (K) 


495.000 

490.000 

485.000 

480.000 

475.000 

470.000 

465.000 

460.000 

455.000 

450.000 

445.000 

440.000 

435.000 

430.000 

425.000 

420.000 

415.000 

410.000 

405.000 

400.000 

395.000 

390.000 

385.000 

380.000 

375.000 

370.000 

365.000 

360.000 

355.000 

350.000 

345.000 

340.000 

335.000 

330.000 

325.000 

320.000 

315.000 

310.000 

305.000 

300.000 

295.000 

290.000 

285.000 

280.000 

275.000 

270.000 

265.000 

260.000 

255.000 

250.000 

245.000 

240.000 

235.000 


Theta- 1 
Gamma- 1/2 

675.848 

230.877 

142.184 

104.542 

84.037 

71.460 

63.333 

56.368 

50.779 
46.332 
42.690 
39.639 
37.036 

34.779 
32.796 
31.035 
29.454 
28.022 
26.71-6 
25.516 
24.405 
23.373 
22.407 
21.501 
20.646 
19.836 
19.067 


Table  IV.  Average  Lamellar  Thickness  (£)  vs.  Temperature  (K) 

Gamma 

TEMP.  (K)  0.5  //  60  0.5  //  100  0.5  II  ISO 


485.000 

235.303 

235.303 

235.303 

480.000 

179.781 

179.781 

179.781 

475.000 

146.556 

146.556 

146.556 

470.000 

124.507 

124.507 

124.507 

465.000 

108.867 

108.867 

108.867 

460.000 

97.253 

97.253 

97.253 

455.000 

88.342 

88.342 

88.342 

450.000 

81.344 

81.344 

81.344 

445.000 

75.762 

75.762 

75.762 

440.000 

71.267 

71.267 

71.267 

435.000 

tt7.641 

67.641 

67.641 

430.000 

64.735 

64.735 

64.735 

425.000 

62.454 

62.454 

62.454 

420.000 

60.723 

60.743 

60.743 

415.000 

59.214 

59.577 

59.584 

410.000 

57.306 

58.874 

59.005 

405.000 

54.856 

58.337 

5B.984 

400.000 

52.149 

57.582 

59.533 

395.000 

49.469 

56.411 

60.296 

390.000 

46.971 

54.852 

60.919 

385.000 

44.708 

53.035 

61.120 

380.000 

42.683 

51.095 

60.800 

375.000 

40.874 

49.136 

60.003 

370.000 

39.252 

47.220 

.  58.842 

365.000 

37.791 

45.385 

57.434 

360.000 

36.466 

43.648 

55.882 

355.000 

35.253 

42.015 

54.261 

350.000 

34.136 

40.486 

52.625 

345.000 

33.100 

39.056 

51.009 

340.000 

32.131 

37.719 

49.436 

335.000 

31.219 

36.468 

47.918 

330.000 

30.358 

35.295 

46.462 

325.000 

29.541 

34.194 

45.072 

320.000 

28.766 

33.159 

43.746 

315.000 

28.028 

32.183 

42.485 

310.000 

27.324 

31.262 

41.286 

305.000 

26.652 

30.390 

40.145 

300.000 

26.009 

29.564 

39.059 

295.000 

25.392 

28.778 

38.025 

290.000 

24.799 

28.031 

37.040 

285.000 

24.229 

27.318 

36.101 

280.000 

23.6B1 

26.637 

35.204 

275.000 

23.152 

25.985 

34.346 

270.000 

22.641 

25.360 

33.526 

265.000 

22.147 

24.760 

32.740 

260.000 

21.669 

24.184 

31.986 

255.000 

21.206 

23.628 

31.263 

250.000 

20.756 

23.093 

30.567 

245.000 

20.320 

22.576 

29.898 

240.000 

19.896 

22.076 

29.253 

235.000 

19.484 

21.593 

28.631 

Table  V.  Average  Lamellar  Thickness  (8‘ 
vs.  Tempera ture  (K) 

TEMP.  (K)  Gamma-0.90 


485.000 

236.013 

480.000 

180.939 

475.000 

148.300 

470.000 

127.027 

465.000 

112.435 

460.000 

102.279 

455.000 

95.475 

450.000 

91.672 

445.000 

91.275 

440.000 

96.119 

435.000 

112.616 

430.000 

117.625 

425.000 

109.730 

420.000 

103.882 

415.000 

99.316 

410.000 

95.563 

405.000 

92.353 

400.000 

89.529 

305.000 

86.992 

390.000 

84.676 

385.000 

82.538  ' 

380.000 

80.545 

375.000 

78.673 

370.000 

76.905 

365.000 

75.225 

360.000 

73.622 

355.000 

72.087 

350.000 

70.612 

345.000 

69.190 

340.000 

67.817 

335.000 

66.486 

330.000 

85.194 

325.000 

63.938 

320.000 

62.714 

315.000 

61.519 

310.000 

60.352 

305.000 

59.210 

300.000 

58.090 

295.000 

56.992 

290.000 

55.913 

285.000 

54.853 

280.000 

53.809 

275.000 

52.782 

270.000 

51.769 

265.000 

50.770 

260.000 

49.784 

255.000 

48.810 

250.000 

47.847 

245.000 

46.895 

240.000 

45.953 

235.000 

45.021 

Table  VI.  Average  Lamellar  Thickness  {%)  vs.  Temperature  (E) 
TEMP.  (K)  Psi-0.90 


485.000 

237.166 

480.000 

182.177 

475.000 

149.552 

470.000 

128.225 

465.000 

113.507 

460.000 

103.129 

455.000 

95.962 

450.000 

91.560 

445.000 

90.139 

440.000 

93.098 

435.000 

105.777 

430.000 

160.924 

425.000 

OO 

